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Abstract 

Let v be a bounded function with bounded support in R d , d > 3. 
Let !,!/£ R d . Let Z(t) denote the path integral of v along the path 
of a Brownian bridge in R d which runs for time t, starting at x and 
ending at y. As t — > oo, it is perhaps evident that the distribution of 
Z(t) converges weakly to that of the sum of the integrals of v along the 
paths of two independent Brownian motions, starting at x and y and 
running forever. Here we prove a stronger result, namely convergence 
of the corresponding moment generating functions and of moments. 
This result is needed for applications in physics. 
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1 Introduction 



Suppose that (X s , s > 0) is a stochastic process taking values in R d , and v 
is a Real-valued function on R d . Define the path integral Z(t) by 



If v is the indicator function of a set Z(t) is the occupation time of that set 
by (X s ) up to time t. 

If (X s ) is a suitable time- homogeneous Markov process, such as a Feller 
process, then the distribution of Z(t) can be studied via the Feynman-Kac 
formula. See for example Williams [8]. For example, if (X s ) is Brownian 
motion and v is the indicator function of the unit ball, the distribution of 
Zoo was obtained by Cieselski and Taylor [3] (see also [8]). In the Markov 
process literature the path integral Z(t) is known as an additive functional. 

We study here the path integral Z(t) given by (1), where the process 
(X s , < s < t) is a Brownian bridge, and therefore is not time-homogeneous. 
Loosely speaking, a Brownian bridge (X s ,0 < s < t) is a Brownian motion 
conditioned to take some fixed value y at time t. 

Our motivation for this work came from the study of the heat equation 
with added potential v(-); that is the equation 



(where the Laplacian y 2 acts on functions of x at fixed t). This equation 
is closely related to the Schrodinger equation. In quantum Statistical Me- 
chanics, in [6] and [7] for example, one studies the second virial coefficient 
for a gas consisting of particles with a two-body interaction potential v, 
via the fundamental solution G(t;x,y) = G y (t;x) to the associated Bloch 
equation (2); that is, the solution / = G y to (2) with the initial condition 



By the Feynman-Kac formula, one can show G(t; x, y) is equal to the 
Brownian transition density (2irt)~ d / 2 e~\ y ~ x \ K 2t \ multiplied by the moment 
generating function (m.g.f.) of a random variable of the form Z(t) given by 
(1) where (X s , < s < t) is a Brownian bridge running from x to y in time 
t (see [7], equation (29) for details). It is this m.g.f. which we study here. 

Another motivation arises in the simulation of molecular dynamics for a 
chemical reaction. See Clifford et al. [2]. Suppose two particles (molecules) 




(1) 



df/dt = [(1/2) v 2 -v(x)]f(t; x) (t > 0, x e R d ) 



(2) 



G y (0;x) = 5(x - y). 
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execute independent Brownian motions, so that one particle executes a Brown- 
ian motion (W s ) relative to the other. Suppose the particles react (for exam- 
ple by mutual annihilation) at a rate given by a function v of their separation 
vector W s ] that is, if r denotes the reaction time, 

PW <t]= Eexp{- f v(W s )dx}. (3) 
Jo 

In a simulation of the system of two particles (maybe as part of a much 
larger system), it is natural to simulate the value of W s , say at integer values 
of s, and to interpolate between successive integer times n and n + 1 by 
taking the path of (W s ) between those times to be a Brownian bridge with 
endpoints given by the simulated values of W n and W n+ i. Given that no 
reaction takes place before time n, and given the simulated values of W n and 
W n+ i, the probability of reaction before time n + 1 is again the m.g.f. of a 
path integral of v along the path of a Brownian bridge. 

This approach was used by Clifford and Green [1], although their reaction 
mechanism is different from the one discussed here. The Brownian bridge 
could arise similarly in simulations of an integrated cost function in Control 
theory. 

In the physical examples discussed above, the process (X s ) of (1) is typ- 
ically a Brownian bridge in three dimensions. In three or more dimensions, 
Brownian motion is transient, and if v has bounded support, denoted A say, 
then for large t one might expect the Brownian motion to spend most of its 
time outside A, leaving A shortly after time and returning shortly before 
time t (where 'shortly' is relative to t). Since the time- reversed process is 
also a Brownian bridge, for large t one might expect the distribution of Z(t) 
to approximate to that of 

roc roc 

J o v(W x (s))ds + j o v(W^s))ds, (4) 

where W x and W are independent unconditioned Brownian motions starting 
at x and y (the endpoints of the bridge). Theorem 1 below makes this intu- 
ition precise; we go beyond convergence in distribution, proving convergence 
of m.g.f. s and of moments. 

2 Statement of results 

Assumptions. We assume throughout that d > 3 and that v is a bounded, 
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Borel measurable function on R d with bounded support. Let x(t) G R d and 
y(t) G R d be defined for t > 0. For each t > let the random variable 
be given by the path integral (1) of v, along the path (X.,0 < s < t) 
of a ci-dimensional Brownian bridge running for time t with X = x(t) and 
X t = y(t). We view {Z(t),t > 0} not as a stochastic process, merely as a 
family of probablity distributions indexed by t. 

For x G let (W x (s), s > 0) be a Wiener process (unconstrained Brown- 
ian motion) with Wc(0) = x. For t G [0, oo], define 

Y x (t)= fv{W x {s))ds. 
Jo 

Write Y x as an abbreviation for Y x (oo). 

Before stating our main results on Brownian bridge path integrals, we 
state some facts about the Brownian motion path integral X x . 

Lemma 1 There exists a > such that for \a\ < ao, 

Eexp{a \v(W x (s))\ds} < oo, x G R d . (5) 
Jo 

Lemma 2 For every real a, Ee aYx is either finite for all x or infinite for all 
x. In the former case, Ee aYx is a locally bounded, continuous function of x. 
Also, for all nonnegative integers k, E(Y x ) k is continuous in x. 

Lemma 3 Suppose v is nonnegative and not almost everywhere zero. Then 
there exists a.\ G (0, oo) such that 

Ee aYx < oo, a < a±, x G R d ; (6) 
Ee aY * = oo, a > cci, x G R d . (7) 

We can now state our main results. All limits below are taken as t — > oo. 
This limiting procedure is motivated by problems arising in [7]. Using the 
scaling property of Brownian motion, one can re-state the results in terms of 
a limiting regime where the time for which the Brownian bridge runs remains 
fixed, and the range (support) of v shrinks. 

Convergence in distribution is implied by the statements (8) and (9) below 
on convergence of m.g.f.s on a neighborhood of (see [5]). Our results on 
convergence of m.g.f.s are stronger than convergence in distribution, and are 
used in [7]. 
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Theorem 1 Suppose x(t) — > x G -R d , and — > y G Tnen eac/i 

moment of Z{t) converges to the corresponding moment of Y x + Yy, where 
Yy is defined to be a random variable, independent of Y x , with the same 
distribution as Y y . Also, for a given by Lemma 1, 

Ee aZ{t) -> Ee^'+W = Ee aY *Ee aY \ a G (-a , a ). (8) 

// v is nonnegative, then the convergence in (8) holds for every a G (— oo, ai), 
where a.\ is given by Lemma 3. 

Theorem 2 Suppose x(t) — > x G -R d , w/w/e — *■ °°; an d \y(t)\ 2 /t re- 
mains bounded. Then every moment of Z(t) converges to the corresponding 
moment ofY x , and for \a\ < ao, 

Ee aZ ® -> Ee aY ", a G (-a ,a ). (9) 

If v is nonnegative, then (9) holds for every a G (— oo,ai). 

To prove the Theorems, we need the following extension of the continuity 
results of Lemma 3. 

Lemma 4 (a) Suppose (x n ) and (t n ) are sequences in R d and (0, oo) respec- 
tively, with x n — * x G R d and t n — > oo as n — > oo. Then each moment of 
Yx n {tn) converges to the corresponding moment ofY x . Also, 

E exp(aY Xn (t n )) -> E exp(aY x ), a G (-a ,a ). (10) 

If v > on R d , then (10) holds for alia G (— oo, cti). 

(b) Suppose now that \x n \ — > oo. Tnen Ee aYxn( - tn ^ — > as n — > oo, /or 
any /raed a G (— a , ao) and /or any /raed a G (— oo, ai) if v is nonnegative. 

3 Proof of Lemmas 

Notation. Let q(t; x) (t > 0, x G -R d ) denote the transition density of d- 
dimensional Brownian motion; that is, q(t;x) = (2nt)~ d / 2 exp(— \x\ 2 /(2t)). 
For s = (si,...,s fc ) G (0, oo) fc and z = (zi,...,z k ) G (i? d ) fc , let q%(s, z) 
denote the joint probability density function at (zi, . . . , Zk) for the positions 
at times (s 1? . . . , s fe ) of a d-dimensional Brownian motion starting at x. Let 
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qtk(s, z) denote the corresponding joint density of a (/-dimensional Brownian 
bridge (X s , < s < t) starting at x and finishing at y. 

For t G (0, oo], write T(k; t) for the 'triangular' set {s = (s l5 . . . ,s k ) : < 
Si < • • • < s k < t}. For s G T(k; oo), we have 

k— 1 

gJ(s,z) = n?(A^;A^). (11) 

where we set so = 0, = AjS = Sj+i — Sj and A^z = Zj+i — Zj (0 < j < 
k — 1). Similarly for the Brownian bridge; for s G T(/c; t), 

«& W (s, *) = [f[ g(A jS ; A,-*)]/?*!/ - x; t), (12) 

j=0 

where we set s k +i = t, and x^+i = y. 

For z = (z!, z 2 , • • • , z k ) G (i? d ) fc , set t> fc (z) = v(z 1 )v(z 2 ) • • ■ v(z k ). 

Proof of Lemma 1. We have 

E {£ H W I{S ))\ is)" IH = / T<to) ( W z)|,*(s, .)*) *. (13) 
Since is bounded and has bounded support, there exists finite K\ for which 

/'OO /■ 

/ / |v(z)|?(s;z-y)dzds < Xi, |/6^. (14) 

This implies that the expression in (13) is at most K\\ hence for \a\ < 
Ki 1 , the expansion of the m.g.f. (5) as a series of moments is absoloutely 
convergent. 

Proof of Lemma 2. Set f(x) = Ee aYx . Suppose for some x G R d that 
f(x) < oo. We must show that f(y) < oo for all y G R d . 

Let B be a ball of radius r, centered at x. Take s > 0, and for all y <E B, 
set Ty = inf{t > : W y (t) ^ B}, and = min(r y , s). Define the measure \i y 
on 95, by /^(F) = P[W y {a y ) G £]. Let K = sup{|u(;s)| : z G R d }. By using 
the Markov property for W x at time s, 

f(x)>e~ Ks [ f(z)q(s;z-x)dz (15) 
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and by the strong Markov property at time a x , 

f(x)>e~ Ks [ f(z)pi x (dz). (16) 

JdB 

Suppose \y — x\ < r/2. By the strong Markov property at time a y , 

f(y) < e Ks (J b f(z)q(s; z - y)dz + f(z)/, y (dz^j . (17) 

Since for s fixed, q(s; z — y)/q(s; z — x) is uniformly bounded on z G B, and 
\i y has bounded Radon-Nikodym derivative with respect to fi x , it follows that 
f(y) < oo, and indeed / is bounded on {y : \y — x\ < r/2}. Since r was 
arbitrary, / is everywhere finite and is locally bounded. 

We now show / is continuous. In the above, let the ball B have radius 1. 
Let M be the supremum of / over dB, which is finite. The second integral 
in (17) is at most MP[r y < s], so if y n — > x, then by (17), 

limsup f(y n ) < e Ks f(z)q(s; z - x)dz + MP[r x < s]) . (18) 

Now combine (18) with (15), and make s — > 0, to obtain limsup f(y n ) < f(x). 
A similar argument gives us liminf f(y n ) > f(x). 

To show E(Y x ) k is continuous in x for nonnegative integer k, use induction 
on k. Set 

f k (x) = {Ee aY * - X: a*E(Y x y/j\)/a k+1 . 

j=0 

Then f k is continuous in x by the inductive hypothesis, and as a — > 0, fk(x) 
converges to E(Y x ) k+l / (k + 1)!, uniformly in x (see the proof of Lemma 1). 

Proof of Lemma 3. If v is nonnegative, then Ee aYx is clearly monotone in 
a, and Lemmas 1 and 2 imply the existence of a± G (0, oo] satisfying (6) and 
(7). Also, a± is finite; indeed, a straightforward geometric series argument 
[4] shows that for given x, Ee aYx is infinite for large enough a. Since we shall 
not use the fact that cci < oo here, we omit details. 

Proof of Lemma 4- (a) First suppose v is nonnegative. By Lemma 2, 

E exp{aY Xn (t n )} < E exp{aY x J ^Eexp{aY x }. (19) 

Also, for fixed to < oo, the series (E\aY Xn (t )\ k /k\,k > 0) is dominated (for 
all n) by the summable series ((aKt ) k /k\, k > 0). Each moment of Y Xn (t ) 
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converges to the corresponding moment of Y x (to); to see this, write each 
moment as a multiple integral with respect to the transition density q^ n or 
q%, and use a domination argument. Hence, 

Eexp{aY Xn (t )} — > E exp{aY x (t )} , (20) 

and by making t — > oo, 

liminf i?exp{o;y En (t n )} > explain} 

so that by comparison with (19) we have (10). 

Now drop the assumption that v is nonnegative. If |a| < cko, then (10) 
would hold if we replaced v(-) by \v(-)\ in the definitions of Y Xn (t n ) and Y x . 
Hence by obvious term-by-term estimation on the difference of the power 
series, \Ee aYx ^ — Ee aYx \ — > locally uniformly in x, and (10) follows. The 
proof of convergence of moments is similar. 

(b) This is a routine application of the strong Markov property, using the 
fact that v has bounded support and Ee aYx is locally bounded in x. 

4 Proof of Theorems 

Take an arbitrary increasing Real-valued function u(t),t > with u(0) = 0, 
and with u(t) — > oo and u(t)/t — > as t — > oo. The idea of this proof is to 
compare the moments of Z(t) with those of Y x ^(u(t)) + YL t Au{t)). 

Proof of Theorem 1. It suffices to prove the results along an arbitrary 
sequence (t n ) with t n — > oo. Write x n for x(t n ), y n for y(t n ), and w n for 
u(t n ). Suppose first that v is nonnegative. 

Let Ti(k; t n ) denote the set {s e T(fc; t n ) : Sj G (0, w n ) U (t n — u n , oo), 1 < 
j < fc}, and let T 2 (/c; t n ) denote the set T(k; t n ) \ T\{k\ t n ). Then 

E(Z(t n ) f/k\ = f f CT(s; z) Vfc (z)dzds (21) 

JT(k;t n ) J (R ) 

= h + h (22) 

where 7j denotes the integral over Ti(k;t n ) (i = 1, 2). 

Now consider the sum of independent unconstrained Brownian path in- 
tegrals Y Xn (t n ) + Yy n (t n ). By the Binomial theorem, 

E(Y Xn (t n )+Yi( tn )) k /k\ = i:{E[Y x stn)\ j /jm[Y yn {tn)\ k - j i{k-m 

3=0 
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x / v j (z)q^ n (s;z)dzds. (23) 

Indeed, the expression inside the braces is less than E(Y yn {t n )) k ~i /{k — 
since s' has been integrated over T(k — j; t n — sj) (sj is the last component 
of s), rather than T(k — j; t n ). 
By the substitutions 

Zj +r = ^fc+l-(j+r)? S j+r — tn — s k+l-(j+r) ( r = l,2,...,fc j), 

the right hand side of (23) equals 

k 

S Ln ^ /D ^ k V ki^)q £ j n {.S U ...,S j ]Z 1 ,...,Z j )x 
^ Q JT(k;t n ) J(R d ) k 

Xql-j(tn — s k, tn — Sfc-l, . . . ,t n — Sj+i, Zfc, 2fc-i, • • • , Zj +1 )dzds. (24) 

On inspection, we find that by (12), (11) and the symmetry of q(t;x) in x, 
the integrand in (24) equals 

u fc( z )C",fc"( s i' •••>**; ••• > ^fc)Qnj(s, z) 
where the density ratio Q n j(s, z) is defined by 

QnJs , ,) = ife^ = f^f exp f !^!! - . (25) 



g(Aj-s; Ajz) \ t n J \ Ajs t n 

Hence, 

k 

E(Y Xn (t n )+Y> n (t n )) k /k\ > ]T L , [qlT(s,z)QnAs,z)vk(z)dzd S 

(26) 

A similar calculation for Y Xn (u n ) +Yy n (u n ) gives rise to a similar formula to 
(23), with s' now integrated over T(k—j; u n ) and t n replaced by u n elsewhere, 
and with equality this time, provided n is so big that u n < t n /2. By use of 
the same change of variables as before (that is, running time backwards from 
t n , not u n ), we obtain 

E(Y Xn (u n )+Y' n (u n )) k /k\ = 
9 



k 

= E L fu , L, k ^(s)CT(s, z)Q nJ (s, z)v k (z)dzds, (27) 

j = Q JTl(k;tn) J (-TL ) 

where lj denotes the indicator function of the set {s e 7\(A;; t„) : Sj < «„ < 

t n ~ Un < Sj+l}- 

Take the summation inside the integral in (27). For each s e Ti(k;t n ), 
the indicator function lj(s) is 1 for exactly one value of j. For that value of 
j, and for Vk(z) ^ 0, the ratio Q n j(s, z) is close to 1 when n is large, because 
of the assumptions on x n , y n , and u n . Indeed, we have 

E{Y Xn {u n ) + Yl{u n )f/k\ ~ h asn^oo, (28) 

in the sense that the ratio of the two sides approaches 1, uniformly in k. 

Turning to I 2 , we have that for sufficiently large n, exp(— \y n — x n \ 2 /t n ) > 
1/2, and hence 

EQ nj (s 1 z)>(l/2)(fc + l) 1 - d/2 ) (29) 
3=0 

since by Jensen's inequality, J2j=o(Ajs/t n ) d / 2 is minimised by setting sj = 
t n /(k + 1), for all j. So for large enough n, 

I 2 < 2(k + l)^ 2 - 1 W / q* :: t (s, z)Q n ,,(s, z) Wfc (z)dzds (30) 

j'=0 2 ' >*™) ' ' 

which implies (by (27) and (26)) that 

I 2 < 2{k + l) d / 2 - 1 D{n,k) (31) 

where we set 

D(n, k) = [E(Y Xn (t n ) + Y> n (t n )) k - E(Y Xn (u n ) + Y^(u n )) k ]/k\ (32) 

By Lemma 4, D(n, fc) — > as n — > 00 ; combining this with (28) and (22), 
we obtain the desired convergence of moments. 

Now consider I 2 as a function of k. Suppose < a < a±, and choose 
a 2 G (a, a^. For a suitable choice of k , we have 

OO fco OO 

E W < 2 £ a*(* + l)" 72 " 1 ^, *)) + E ^(D(n, k)). (33) 

fc=0 fc=0 k=k 
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The first sum converges to zero as n — > oo. The second is bounded above by 

E exp{a 2 (Y Xn (t n ) + ^ (*„))} - E exp{a 2 (Y Xri (u n ) + Yy n (u n ))}. 

By Lemma 4, the last expression converges to zero. This, together with the 
uniform convergence in (28), implies the desired convergence (8) of m.g.f.s 
on < a < cki, for v nonnegative. 

Now drop the assumption v is nonnegative. The proof of (28) is un- 
changed. Also, (31) would hold if v were replaced by \v\ in the definitions of 
li and D(n,k); hence I/2I — ► 0, and for a G (— ao,ao), J2k>o \ a \ k \h(k)\ — > 
as above, giving us the desired convergence (8) of m.g.f.s. Finally, this implies 
convergence in distribution, so (8) holds for a < when v is nonnegative. 



Proof of Theorem 2. We use notation from the previous proof. First 
consider the case that \y n \ — * 00, with \y n \ 2 /t n — > 0. Then convergence in 
the desired sense of Z(t n ) to Y x follows by the proof of Theorem 1 (now using 
Lemma 4(b)). 

Now consider the case that \y n \ 2 /t n is bounded away from and from 
00. First assume v is nonnegative. Partition the set T(k;t n ) as T 3 (k;t n ) U 
T4(k;t n ), where we set 

T 3 (k; t n ) = T(k; u n ), T 4 (k; t n ) = T(k; t n ) \ T 3 (k; t n ). 

For i — 3,4, let U denote the integral in (21), restricted to Tj(/c;t„); that is, 

Ii = I f CT( g ; z> k (z)dzds. (34) 

JTi(k;t n ) J(R d ) k 

So E(Z(t n )) k /k\ =I 3 + h. 

A similar (easier) argument to the derivation of (27) gives us 

E(Y Xn (u n )) k /k\ = I f q*:;t(s,z)Q n , k (s,z)v k (z)dzds, (35) 

JT(k;u„) J(R d ) k 

where <3n,/c( s , z ) is given by (25), with j = k and A k z = y n — z k and A k s = 
t n — s k as before. Since \y n — z\ 2 /\y n \ 2 — > 1 as n — > 00, locally uniformly in 
z, we can write 

Q„ lfc (s,z) = (A k s/t n ) d/2 exp{(\y n \ 2 /t n )[(^-)g n (s,z) - h n (s,z)}} 
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where g n and h n converge to 1, uniformly on s G T(k, t n ) and z in the support 
of v k . 

Observe that Q„ jfc (s, z) — > 1, uniformly on s G T(k;u n ) and z in the 
support of Vk- Therefore, ^3 ~ E(Y Xn (u n )) k asn-> 00, uniformly in k. 

To estimate I4, note that \y n \ 2 /t n is assumed to be bounded away from 
and ooasn^ 00, so we can find 5 > and Uq G Z for which Q nt k(s, z) > 5, 
for all s G T(/c; t n ), z G Supp(ffc), n > n and > 0. Thus for n > n , 

/4 < / Cfe"(s; Z)<5n,fe(s, zH(z)dztfe. 

A similar expression to (35) holds with w n replaced by t n , so 

/ 4 < 5- l [E(Y Xn (t n )) k - E(Y Xn (u n )) k ]/k\ 

which implies that as n — > oo, ^4 — > and X)fc>o ® k h(k) — * (0 < ct < 
ax). Combining the estimates obtained for I 3 and 7 4 , we obtain the desired 
convergence results for v nonnegative and a G [0, The full result follows 
as in the proof of Theorem 1. 
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